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Abstract Experimental measurements of the response of is the understanding of the role of fluctuations in the mod- 
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r\^\ • a two dimensional system of plastic beads subjected to ern theory of phase transitions £Q. Similarly, one expects 

o 

^ ■ steady shear are reported. The beads float at the surface nonequilibrium fluctuations to play an equally important 

c3 . 

of a fluid substrate and are subjected to a slow, steady- role in transitions in driven systems, such as the proposed 

shear in a Couette geometry. The flow consists of irregular jamming transition as a function of external stress [21131 

intervals of solid-like, jammed behavior, followed by stress In driven systems, the fluctuations of interest are typi- 

relaxations. We report on statistics that characterize the cally not thermal in nature. Instead, they are the result of 

J> , stress fluctuations as a function of several parameters in- the external driving, and this raises important questions 

C> eluding shear-rate and packing density. Over a range of and presents a number of challenges. As a starting point, 

m . 

' densities between the onset of flow to the onset of buck- it is useful to characterize the statistics of these "ather- 
t-H ' 

' ling (overpacking) of the system, the probability distribu- mal" fluctuations and ask how similar such statistics are 
■ tion for stress fluctuations is essentially independent of the to equivalent measures of thermal fluctuations. Then, one 

S : 

I packing density, particle dispersity, and interaction poten- can consider the impact athermal fluctuations on the sys- 

Sh , tial (varied by changing the substrate). Finally, we com- tern's dynamics. 
O ■ 

^ ' pare the observed stress fluctuations with those observed 
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A general feature of fluctuations in many driven sys- 



! > ' in other complex fluids. . . , r ^ . r . „ 

. terns is the occurrence of non-Gaussian features, typically 



an asymmetric shape with exponential tails, in probability 
distributions of the fluctuating quantity. It is interesting 
to ask under what conditions such distributions are uni- 
versal. This question has been studied in the context of 

For equilibrium systems, it is well understood how fluctua- , . , „ . . n . ,, . . „ . , m ,,, 

turbulent fluid flows for thermally driven fluids 5 (ther- 
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Introduction 



tions are controlled by temperature. Especially important 



mal convection) and electrically driven fluids [HIE] (elec- 
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,,. , , ^ , , ,,. , , r~. . troconvection). More recent work has considered this ques- 
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Two different approaches for characterizing fluctuations in 
the dissipated power were used. The first is based on the 
generalized fluctuation-dissipation theorem of Gallavotti 
and Cohen, also referred to as the Fluctuation-Theorem 
A key aspect of the FT theorem is the occurrence 
of rare negative fluctuations in power. This has been in- 
directly tested in fluid systems and more recently 
direct measurements have been made in a shaken, granu- 
lar system The second approach focuses on the form 
of the probability distribution function, or distribution 
function for short. One prediction is that the generalized 
Fisher-Tippit-Gumbel (gFTG) distribution function de- 
scribes power fluctuations |10| . Experiments using electro- 
convection have demonstrated an interesting evolution of 
the distribution function from Gaussian behavior to be- 
havior that is consistent with the gFTG distribution as 
the distance from equilibrium is varied \7\. 

Another class of systems that have been used heav- 
ily for the study of nonequilibrium fluctuations are soft- 
matter systems. A striking observation that arises from 
these studies is the similarity between the observations of 
fluctuations in foam |11II12II13II14II15) . granular materials 
|16U17lll8ll5) . and colloidal suspensions |T^]. In the exam- 
ples cited here, the fluctuations in stress under conditions 
of steady applied strain rate were studied. Similarities in 
various measures of the fluctuations exist despite qual- 
itative differences in the "microscopic" structure of the 
materials. For example, foam consists of gas bubbles with 
liquid walls. For flowing foam, viscous dissipation in the 
fluid walls is the major source of energy loss. For compar- 
ison, the dissipation in granular systems is usually due to 



dry friction between particles. Another potentially signif- 
icant difference is the the degree to which the particles or 
bubbles are packed. Because of the solid nature of granu- 
lar particles, they are generally studied at packing densi- 
ties near or below random close packed. In contrast, foam 
"melts" at packing densities near random close packed, 
and studies focus on packing densities greater than the 
melting transition. Given these "microscopic" differences, 
it is worth commenting on some of the similarities that 
are observed in previous studies. 

For foam, there is both a rich history of theoretical 

studies |11II12II14II15II13) that deals with the issue of fluc- 
tuations. The focus is on predicting the probability distri- 
bution of energy (or stress) drops, referred to as energy 
(stress) drop distributions. (Typically, simulations focus 
on energy drops because it is relatively easy to compute 
the total energy of the system. In contrast, experiments 
tend to focus on measurements of stress, as this is the more 
accessible quantity.) The focus on the energy (or stress) 
drop distributions was originally motivated by potential 
connections to self-organized criticality and the observa- 
tion of power-law behavior [50]. In this context, the models 
of foam can be divided into two main classes: models that 
predict power-law distributions for stress or energy drops 
[21)112111 22, 23. 24, and models that predict a definite aver- 
age stress drop size [25 , 26 27 28 . Experimental evidence 
from model two-dimensional systems generally supports 
a distribution with a well-defined average |11U12II14II15| . 
However, there is indirect evidence for power-law distribu- 
tions under conditions of quasi-static shear |13) . One issue 



that has not been thoroughly explored in experiments is 
the possibility of the development of power-law behavior 
near the melting transition. There is some evidence in sim- 
ulations exhibiting a well defined average energy drop that 
the average energy drop increases as one approaches the 
melting transition. This might indicate that the distribu- 
tion is approaching a power law |28| . 

For granular materials, there have been limited mea- 
surements of stress drop distributions. Behavior similar to 
that observed in foam is reported for a system of elastic 
disks that are slowly sheared JHj- In general, the distribu- 
tion is consistent with the existence of a characteristic size 
for the stress drops. The system reported on in Ref. [T%| 
exhibits a continuous transition at which stress chains first 
extend across the system Qjj. Near this transition point, 
there is evidence of possible power law behavior for the 
distribution of stress drops. 

More commonly, the focus in granular systems has 
been on the probability distribution for the stress itself 
|16U17ll5] . as was also studied in colloidal systems |19) . 
(For the purposes of this paper, this will be referred to as 
the stress distribution, in contrast to the stress drop dis- 
tribution.) For all of these systems, one surprising feature 
is how large the stress variations can be [161119) . In one 
case, stress was measured normal to the flow in a local- 
ized region of the system. For different particle sizes, small 
variations in a general Maxwellian stress distribution were 
observed, but all of the distributions had essentially expo- 
nential tails for stress greater than the average stress |16| . 
The nature of the stress distributions appeared to depend 
strongly on the nature of the stress chains in the granular 
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material Another case measured stress distributions 
in a colloidal system as a function of strain rates. For low 
strain rates, gaussian distributions in stress were observed. 
For higher strain rates extreme value distribution statis- 
tics were measured with a well defined maximum value 

As a bridge between foam and granular systems, we 
have carried out studies of fluctuations using a modifi- 
cation of the bubble raft [2H]: a bead raft (plastic beads 
floating on the surface of a fluid substrate). The plastic 
beads provide a connection with granular systems. How- 
ever, their individual density is such that they are mostly 
submerged below the fluid substrate. Therefore, there is 
a fluid layer between the beads, analogous in some ways 
to the fluid walls in a foam. In this regard, the system 
is closest to a two-dimensional colloidal system. However, 
in contrast to typical colloidal systems used to study flow 
(see for example, Ref. ^1) the particles used here are 
much larger (mm in size), so thermal fluctuations play no 
role. Also, because the particles are only partially sub- 
merged, the role of the lubrication of the fluid subphase 
is not well understood yet. 

The bead raft is essentially a two dimensional sys- 
tem with free boundary conditions perpendicular to the 
surface. The system has an easily adjustable dispersity, 
packing fraction, strain rate, and interaction potential. We 
report on measurements of stress during steady shear in- 
cluding average stress values (which determine the flow 
regime for a given strain rate), stress distributions, and 
stress drop distributions. By measuring both the stress 
and the stress drop distributions, we are able to make 
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comparisons with the stress drop distributions observed 
in foam and the stress distributions studies in granular 
matter. Using the stress distributions, we can make indi- 
rect comparison with the power dissipated distributions 
observed in turbulent flows. 

2 

Experimental Details 

In our experiments, we use a centrosymmetric trough that 
can measure stress with a freely hanging torsion pendu- 
lum. This apparatus is described in detail in Ref. |30| . 
It consists of two concentric cylinders oriented vertically. 
The outer cylinder consists of 12 individual "fingers" that 
can be expanded and compressed between 6 and 12 cm to 
adjust the packing fraction of the beads. This barrier can 
be rotated in either direction at a constant angular speed 
in the range of 0.0005 rad/sec to 0.4 rad/sec. The inner 
cylinder is a Teflon rotor that "grips" the material by the 
insertion of eight 5 mm aluminum "paddles" . The rotor 
has a diameter of 2.7 cm and is suspended by a torsion 
wire with a torsion constant, k — 3020 dyne cm/rad. In- 
duced voltage in a coil attached to the pendulum was used 
to determine the angle of the pendulum. When the outer 
barrier is rotated we can measure the torque on the inner 
cylinder and the corresponding stress on the inner cylin- 
der. For the rest of the paper, unless otherwise indicated, 
the stress refers to the stress on the inner cylinder. 

The system consists of spherical plastic beads localized 
to the surface of a fluid by their buoyancy. Their density 
is 0.95 g/cm 3 , allowing each sphere to float at the fluid's 
surface. Unless noted, a bidisperse mixture of beads was 



used: 300 0.25 inch beads and 300 0.1875 inch beads. The 
simplest definition of packing fraction was used: the num- 
ber of beads of each size times the maximal cross sectional 
area divided by the area of the system. At the highest 
compressions, this does not account for possible displace- 
ments perpendicular to the fluid surface. The system was 
thoroughly mixed and a disordered two dimensional sys- 
tem was created in the apparatus between the barriers. 
Because some signs of aging on short time scales were 
found, all measurements were performed after 1 hour of 
rotation in the same direction. Over the time scale of our 
experiments, we found no signs of coarsening or other set- 
tling phenomena common with some agitated polydisperse 
granular systems. 

The substrates used ranged from pure water to pure 
glycerine, and included glycerine/water mixtures. The vari- 
ation in substrate effectively changes the interaction be- 
tween the beads because of changes in the properties of 
the fluid layer between particles. For all cases, coupling of 
the beads to the substrate was negligible. This was tested 
by removing the outer layer of beads from contact with 
the outer barrier. In this case, rotation of the outer bar- 
rier yielded zero stress on the inner rotor, and no motion 
of the beads was detected in this case. 

For many of the experiments, the system was video 
taped to monitor the general behavior of the beads. Ex- 
cept for the highest packing fractions (as discussed be- 
low), there was no discernable motion of the beads out 
of the plane of the fluid surface. However, this does not 
completely rule out that some of the release of stress was 
due to slight motions of the beads in this direction. Such 



releases in stress may be important if quantitative com- 
parison with purely two-dimensional models are carried 
out in the future. However, for the main results of this 
paper, which is focused on the degree of universality of 
the fluctuations, the details of the source of stress release 
are not relevant. 

3 

Results 

Figure ^ shows typical time series of the stress on the in- 
ner rotor as the outer barrier is rotated for three different 
packing fractions. The packing fractions were selected to 
illustrate three qualitatively different behaviors. For the 
lowest packing fraction (solid black line, lowest curve), 
the stress is predominately zero. This is due to the fact 
that under constant shear, the particles are driven toward 
the outer boundary. For sufficiently low densities, this re- 
sults in exactly zero stress on the inner boundary, as no 
particles are in contact with it. As the packing fraction 
is increased, there is a density regime for which particles 
intermittently contact the inner rotor. This results in the 
occasional spikes in the stress. 

The second behavior is illustrated by the red line (mid- 
dle curve). This occurs when there are always particles in 
contact with the inner rotor, but the fluctuations occa- 
sionally result in the stress dropping back to zero. Finally, 
at high enough packing fraction (blue curve, highest av- 
erage stress) , the stress is always nonzero during the time 
for which the flow is observed. The rest of the results focus 
on the average stress value, the probability distribution of 
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the stress, and the stress drop distributions, as calculated 
from time series similar to those illustrated in Fig. ^ 

The first question that is addressed is the dependence 
of the average stress on the rotation rate and the density. 
Figure |2 presents the results for the average stress versus 
rotation rate for three different substrates. For all three 
cases, the behavior is consistent with that of a power law 
fluid. A power law fluid is one in which the viscosity, 77, 
as a function of strain rate, 7 is given by: 77 oc A f( m ~ 1 \ 
For the Couette geometry, 7 oc S7, where Q is the outer 
rotation rate. This gives < a >= rjj oc Q m , From this, we 
find that for glycerine m = 0.53 ±0.02. For 70% glycerine 
30% water mixture, m = 0.35 ± 0.02, and for pure water, 
m = 0.20 ±0.01. 

Before discussing these results further, it is useful to 
consider the average stress as a function of packing frac- 
tion, as illustrated in Fig. [3J For this system, there are a 
number of transitions as a function of packing fraction. As 
already discussed with regard to Fig. ^ for small enough 
packing fractions there is no transfer of stress between 
the inner and the outer barriers. The transition from this 
regime to non-zero average stress is somewhat poorly de- 
fined due to long-time transients. Figure |21 only covers 
densities for which the stress did not go to zero over ex- 
perimental time scales of a few hours. The second transi- 
tion is the "buckling" of the bead raft. This is indicated 
by the change in slope at higher packing fraction and is 
confirmed by imaging the bead raft. 

Before the buckling transition occurs, the average stress 
as a function of density is consistent with a diverging 
power law, suggestive of an approach to the "jamming" 
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transition 2 . Therefore, for the range of density of in- 
terest (below buckling and with nonzero average stress), 
the system is below the jamming transition and should be 
in a purely fluid state. This is consistent with the results 
presented in Fig. |2 from which there is no evidence for 
a yield stress. In contrast, the bubble raft is above the 
jamming transition and clearly exhibits a yield stress |14l 

hhi- 

A critical packing fraction p c for the jamming transi- 
tion is determined by fitting the data to 



The fits are given by the solid lines in Fig. 03 It should 
be noticed that this fit focuses on the divergence of the 
stress at higher packing fraction. This is opposite of the 
definition of the critical packing fraction reported in Ref. [T7] 
where the interest was the initial transition to nonzero 
stress. 

Figures [21 and reveal a dependence of two charac- 
terizations of the average properties of the system on the 
choice of substrate: the exponent of the power-law viscos- 
ity (m) and the critical density for jamming (p c ). All of 
the data can be fit with the same exponent of the stress 
divergence (n — 2). The general trend is for m to in- 
crease as the concentration of glycerine is increased, and 
for p c to decrease as parameters are decreased. The values 
of m were already reported. For p c , some examples are 
p c = 0.74 and p c = 0.78 for the pure glycerine and wa- 
ter substrates, respectively. An interesting feature of p c is 
that it appears to be independent of rotation rate, even for 
the glycerine substrate. The dependence of the parameters 



on substrate presumably reflects changes in the interac- 
tions between the beads due to the changing properties 
of the substrate. This is possible because the substrate 
coats a significant fraction of the beads and forms a thin 
film between particles. Therefore, it is reasonable to con- 
sider variations in substrate composition as variations in 
particle interactions. 

Figure 01 displays the probability distribution for the 
stress for four different packing fractions on a water sub- 
strate. As the packing fraction increases, the average stress 
(< a >) and the variance in stress (5<r) increase. Figure 
shows normalized probability distributions, where we have 
plotted the probability of (a— < a >)/Sa for a number of 
different situations. There are two striking features of the 
normalized distributions. First, for all of the parameter 
values for which we made measurements, the normalized 
distributions essentially collapse on a single curve. Sec- 
ond, the curve is asymmetric, with an exponential tail for 
stresses greater than the average stress. 

In order to make comparisons with other work on stress 
distributions, we attempted to fit the data to three differ- 
ent models. The first model is based on the gFTG distri- 
bution function used to describes power fluctuations |10|: 

11(a) = K{e x ~ eX ) a , (2) 

where x = b(a — s) and a — tt/2. 

The parameters in this function are set by the criteria 
of unit area, zero mean, and unit variance. The plot in 
Fig. uses values for the parameters that are expected 
based on Ref. QH|: s = 0.37 and b = -0.94. The only 
difference is that the constant b is negative because our 



distribution is flipped with regard to the asymmetry mea- 
sured in turbulent systems. 

We also fit our plot to distributions of the form 

P(a) = (a - a*fe- (a - a '^/ <a> (3) 

for n = 1 and 2. These are chosen to compare with the 
results from Ref. [IQ an d a Maxwell-Boltzmann distribu- 
tion, respectively. Here, a* was adjusted to get the best fit 
to the data. Not surprisingly, the gFTG distribution and 
distribution used in Ref. |16| best capture the exponential 
tail for high stress and provide the closest fit to the data. 

Figure shows the distribution of stress drops for 
three packing fractions. The general trend is similar to 
that observed for foam. There is an exponential cutoff in 
the stress drop distribution that results in a well-defined 
average stress drop. Unlike the normalized stress distribu- 
tion, the average stress drop increases as one approaches 
the jamming transition from below. However, the aver- 
age stress drop does not diverge as quickly as the average 
stress itself. Therefore, the average stress drop normalized 
by the average stress actually decreases as a function of 
density (see Fig.Q). 

4 

Discussion 

We have presented results for the response of a plastic 
bead raft floating on a fluid substrate to the application 
of constant rate of strain. The system was selected to pro- 
vide a bridge between similar studies in foam and gran- 
ular materials by combining features of both materials, 
as more traditional collodial systems do. The system was 
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studied in a Couette geometry in which the outer cylinder 
was rotated at a constant angular speed and the stress on 
the inner cylinder was measured. A number of properties 
were measured as a function of packing fraction (or den- 
sity) , applied strain rate (or rotation rate) , and interaction 
strength (or substrate composition). The main focus of the 
study was the behavior of the fluctuations in stress. In this 
regard, two quantities were characterized: the distribution 
of stress and the distribution of stress drops. 

The results for the stress drop distributions are very 
similar to those found in foam |14IH5| and granular ma- 
terial |18| . There is evidence for power-law behavior for 
small stress drops, but there is an exponential cutoff. The 
net result is a well-defined average stress drop. The av- 
erage stress drop does increase as the packing fraction is 
increased, but not as fast as the average stress. There was 
no strong evidence for an approach to a power-law behav- 
ior at either the high or low critical densities. However, this 
might have been due to difficulties with approaching close 
enough to these transitions. For low packing fractions, the 
system became highly inhomogeneous, with beads concen- 
trating toward the outer cylinder. At the high end, there 
was the cross-over to stress releases due to bead motion 
into the third dimension. Future work will be needed to 
explore the detailed behavior near these transitions. 

With regard to the distribution of stress itself, some 
interesting "universal" features are observed. The most 
obvious feature is the exponential tail in the stress distri- 
bution that is consistent with observations in both gran- 
ular systems \fi^ and highly driven fluids [7] . As with the 
driven fluid system, it is intriguing to note that the gFTG 
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distribution function, proposed in Ref. 10 , is consistent 
with the data. The one difference is that the exponen- 
tial tail in our data is for stress greater than the average 
stress. This is similar to the granular systems, in which 
stress is also the measured quantity |T5]. This difference 
can be understood in terms of the connection between 
stress transmitted to the inner rotor and the power dis- 
sipated by the system. Because a constant power is input 
from the outer rotor, when the stress on the inner cylin- 
der is higher than average, the power dissipated by the 
system is lower than average. Additional theoretical work 
is needed to fully understand why, and under what condi- 
tions, the gFTG distribution applies to this system. Also, 
it should be noted that there is still sufficient scatter in 
the data that it is not possible to rule out other distribu- 
tions with exponential tails, such as the function used in 

Ref. nnj. 

One feature of the stress distributions that we were 
not able to access with our system was the extremely high 
strain rate limit. For colloids in this limit ^5], the stress 
distribution developed an interesting structure that in- 
cluded a plateau. This behavior was never observed in our 
system. Future work will attempt to reach the necessary 
high strain rates, as the beads at the air-water interface 
is essentially a colloidal system. 

Finally, the average stress was found to diverge as one 
approached a critical density as the density was increased 
at constant strain rate. This represents a diverging vis- 
cosity for the system. One way to understand this be- 
havior is within the context of the proposed "jamming" 
transition [5100] . Jamming is the crowding of constituent 



particles, disabling their kinetics and further exploration 
of phase space 2,3,4 . When a system develops a yield 
stress or extremely long relaxation time, it jams |31| . The 
jamming phase diagram proposes that one can approach 
the jammed state by either varying density, temperature, 
or shear stress An interesting feature of the jamming 
transition is that for finite size systems in the absence of 
shear, there is a distribution of critical densities that de- 
pend on the particle configurations |32II31II33| . An open 
question is how this distribution of jamming transitions 
impacts the behavior of the system under steady shear. 
It is easy to imagine that the shear causes the system 
to explore configuration space, randomly moving between 
jammed and unjammed states. This exploration of phase 
space would be a source of stress fluctuations. A future 
challenge is to understand the connection between the 
average stress, the fluctuations, and the distribution of 
jammed states. Future experiments exploring the impact 
of system size will be critical to answering this intriguing 
question. 

This work was supported by Department of Energy grant DE- 
FG02-03ED46071, the Research Corporation and Alfred P. Sloan 
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Fig. 1. A plot of stress on the inner cylinder versus time 
during steady rotation of the outer barrier. The characteristic 
stress relaxations are evident. The packing fractions for in- 
creasing average stress are 0.705 (black), 0.74 (red) and 0.78 
(blue). 



Fig. 2. A plot showing average stress values versus rotation 
rate for three different substrates: pure glycerine (•); 70 % 
glycerine/ 30 % water mixture (a); and pure water (■). The 
solid lines indicate the consistency between the data and power 
law behavior. 
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Fig. 3. A plot of the average stress value versus packing frac- 
tion for three different situations: rotation rate of 0.1 rad/s 
with pure glycerine substrate (a); rotation rate of 0.02 rad/s 
with pure glycerine substrate (■); and rotation rate of 0.1 
rad/s with pure water substrate (♦). The solid lines are a fits 
to Eq. Q and they illustrate the diverging nature of the stress 
as the material is compressed. 



Fig. 4. A plot of the probability of a particular stress value for 
beads on pure water with a rotation rate of 0.1 rad/s. The four 
curves correspond to packing fractions of (in order of increasing 
mean): 0.71 (black); 0.74 (red); 0.75 (green); and 0.78 (blue). 
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Fig. 5. The probability of [a— < cr >)/Sa is plotted, where 
8a is the standard deviation of the stress distribution. This is 
referred to as the normalized stress distribution. The different 
substrates and packing fractions are as follows: water substrate 
and packing fraction 0.71 (blue ▼); water substrate and pack- 
ing fraction 0.77 (green A); packing fraction 0.68 and glycerine 
substrate (red •); and a monodisperse system on water susb- 
trate (black ■). The lines are fits to three different possible 
distributions, as described in the text: the gFTG (solid line), 
two-dimensional Maxwell Boltzmann (dashed line), and one- 
dimensional Maxwell Boltzmann (dotted line) 



Fig. 6. A plot of the probability distribution for stress drops 
normalized by the average stress (Aa / < a > for three different 
densities: p = 0.735 (A); p = 0.75 (a); and p = 0.765 (■). 
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Fig. 7. A plot of the average stress drop, normalized by the 
average stress, (< Aa > / < a > versus packing fraction. This 
is for the bidisperse system on water substrate. 



